Dynamics of Multidimensional Secession 
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We explore a generalized Seceder Model with variable size selection groups and higher dimensional 
genotypes, uncovering its well-defined mean-field limiting behavior. Mapping to a discrete, deter- 
ministic version, we pin down the upper critical size of the multiplet selection group, characterize 
all relevant dynamically stable fixed points, and provide a complete analytical description of its 
self-similar hierarchy of multiple branch solutions. 



Dynamical phenomena in which an initially homoge- 
neous population of weakly interacting individual agents 
can disperse, aggregate and form clusters arises in many 
different physical, biological, and sociological contexts. 
Condensation and droplet formation is, of course, a 
well-known example in physics; galaxy formation and 
clustering another. In traffic patterns ||], the real- 
space jams that plague highway driving are, for some, a 
daily reminder of such intrinsic tendencies in correlated 
systems far-from-equilibrium. The formation of swarms 
and herds in zoology [||, or the flocking of birds |,|, 
provide additional illustrations. In these cases, partic- 
ularly, joining the group yields advantages over stand- 
ing out alone, be it by better exploration of food re- 
sources, protection from predators, or easing the aerody- 
namic flow in flight. Nevertheless, sometimes, as in fash- 
ion trends and similar social (or even financial) settings, 
standing apart from the crowd can also be a seed for the 
formation of new groups, splitting off the mainstream, 
though maybe becoming the mainstream themselves later 
on. In these instances, steady-state multiple groups can 
be the norm. Such matters are manifest in recent, though 
now classic implementations of Arthur's variant of the 
El Parol Bar problem 0, as for example, discussed by 
Zhang and Challet ||] , where a multitude of competing 
agents, armed with limited memory strategies, compete 
via statistical Sisyphian dynamics to be in the minor- 
ity group. Interestingly, with stochasticity introduced to 
the decision-making process, Johnson and coworkers |^] 
uncovered a tendency towards self-organized segregation 
within such evolutionary minority games. Subsequently, 
Hod & Nakar |l^] discovered a dynamical phase transi- 
tion in this setting, between 2-group segregation and sin- 
gle group clustering, driven by the economic cost-benefit 
ratio defined in the model. In biological systems, clus- 
tering can appear on multiple scales ||T^, with aggre- 
gation a consequence of dynamic correlations, whether 
they be hidden or explicit. Even so, the complexity 
that arises, for example, in statistical models of evolu- 
tion ijl^, resulting in the formation of species is not, 
per se, self-evident via the direct interplay of mutation 
and selection- the system can dynamically bring itself to 
a critical state. The Seceder Model was introduced, 
initially, to demonstrate that an interative mechanism 
favoring individuality cannot only create distinct groups, 
but also yields a rich diversity of cluster-forming dynam- 
ics. The essential tack was to give a small advantage to 



individuals that distinguish themselves from others. This 
is not unnatural, since in epidemics, for example, genetic 
differences can enhance long-term survival probabilities. 
Likewise, for players in a minority game, distinctness may 
be the advantageous property |f7|,p|js|jl0[| . In this Letter, 
we consider the Seceder Model in its broadest sense. Its 
description is simple enough: within a population of in- 
dividuals, each described by a genotype variable, choose 
a subset from the population and calculate its average. 
From within this selection multiplet, the individual most 
distant from the mean is the parent to be. Create an off- 
spring by taking this parent's value plus a small uniform 
deviate. Finally, replace a randomly chosen member of 
the population by this new offspring. The process is then 
iterated through many generational time-steps. 

Despite the complexity of its segregative dynamics, the 
Seceder Model may be amenable to traditional method- 
ological approaches. For example, the nontrivial scaling 
exhibited by the Seceder envelope, as well as the interest- 
ing time evolution of the group number could be thought 
of as fluctuation-dominated non-classical behavior. In 
this spirit, it is natural to consider the dimensionality 
dependences inherent in the model, with the expectation 
of flnding a simpler, mean-fleld or classical nonequilib- 
rium dynamics within some sector of this larger parame- 
ter space. Some wisdom in this regard may be had from 
the Bak-Sneppen model of punctuated evolution Q, 
wherein a similar, innocuously trim update algorithm 
engenders an extraordinarily rich spatiotemporal dynam- 
ics. Even so, a mean-field limit of this model was sub- 
sequently engineered ||l^,|l^ and further explored 
simplified scaling retrieved by introducing system-wide 
correlations to the interactions. Here, for the Seceder 
Model, we're motivated by similar goals. Clearly, the 
number of parameters is restricted- we have the popula- 
tion size N, which is understood to diverge in the thermo- 
dynamic limit. There is also the size m of the multiplet 
selection group; finally, the dimensionality d of the geno- 
type variable, which determines the nature of the base- 
space through which the population groups mark their 
trajectories, provides an additional degree of freedom. 

We begin by enlarging the selection group from which 
the parent is chosen. Naively, we'd expect the limit 
m N, which introduces increasing cross-correlation 
within the society, to elicit eventually, a mean-field type 
of behavior, if only in the extreme case when m = N, 
when we're averaging over the entire population using 



the societal mean to determine the most distant, repro- 
duced individual. Indeed, this is the case. The surprise, 
however, comes with the abruptness of the transition. 
There is, already, a marked change of behavior as we 
switch from a triplet (m = 3) to a quartet (m = 4) selec- 
tion group. In Figure 1, we show single runs of the d ~ 1 
Seceder Model using multiplet groups m = 3 — 8, within 
an essentially infinite population, N — 512. For to = 3, 
we have trademark Seceder demeanor, with self-similar 
branching characterized by three dominant, but fluctu- 
ating arms, centered about the origin, with ample small- 
scale stochastic structure associated with the transient 
appearance of variously short-lived subbranches. Rather 
than a gradual transition, we find for m — A that the typi- 
cal stable configuration suddenly involves two groups, not 
three. In addition, these two branches exhibit only the 
most modest sorts of fluctuations, as is evident from the 
figure. Increasing the selection group to to = 5 further 
diminishes the fluctuations, but hardly affects the tilt of 
what seems to be the nearly linear divergence of the two 
groups. Next, for to — 6&7, there is, strangely, a dis- 
crete jump to an altogether different, but closely allied 
pair of trajectories. With m = 8(&9), another jump, and 
so it goes with each successive even-odd pair of multiplet 
selection groups. As to A^, the trajectories form an 
extremely well-defined V-shaped wedge, with little fluc- 
tuation at all, as we'd expect of a mean-field limit. Thus 
we see that the dominant dynamic of secession involves, 
for TO > 4, segregation into two evenly populated oppos- 
ing groups with a free interchange of individuals over the 
course of time. Similar self-organized segregation was re- 
ported recently by Johnson et al., |||,|l^ within the con- 
text of an evolutionary minority game. Ensemble av- 
eraging over many realizations, we have systematically 
studied the growth of the population diameter over time. 
Only for triplet selection, to = 3, do we find a fractional 
power-law dependence, the diameter asymptotically scal- 
ing with an exponent very close to 3/4, our measured 
value being 0.74±0.01 for this one dimensional case. 

With the Seceder Model defined as above, the genotype 
space is a continuum. Clearly, discretizing the model al- 
ters no essential features. Indeed, much can be gleaned 
by considering this discrete Seceder Model in its deter- 
ministic limit, wherein the most distinct individual is re- 
produced exactly, rather than yielding a merely approx- 
imate next of kin. One is lead to a set of nonlinear cou- 
pled ODEs , first-order rate equations for the concen- 
tration simplex {xi,X2, ■■■,xb), describing the evolution 
of the discrete set of B genotypes possible within the 
population: xj = Ylfi,...,i^=i(^iii2...i,„^ii^i2---^i,r. - 
for i G (1, B). These equations transform the Seceder 
Model into an evolving chemical reaction system whose 
dynamics are dictated by the law of mass kinetics and the 
constraints of unit dilution flux, possessing some features 
reminiscent of earlier efforts on the hypercycle model, 
and generalized replicator equations p8| . Here, one is 
looking at the stability of an _B— branch solution gener- 
ated by TO— multiplet selection group dynamics. The co- 
efficients are zero unless the genotype/individual is the 
distant outlier- either in isolation, in which case a = 1, 
or as happens occasionally, sharing that distinction with 



009 OOt- 00£ OOe 001 




FIG. 1. Space-time plots of the d = 1 stochastic Seceder 
Model with selection group sizes m = 3 — 8. For m = 3, multi- 
ple groups occur; however, for m > nic — 4, the segregational 
dynamics yields just two repelling clusters, characterized by 
increasing homogeneity in the societal Seceder limit. 

p other individuals, with a ~ l/p. As a practical 
matter, the coefficients, combinatoric in origin, can by 
generated systematically via a multinomial expansion 
{xi + X2 + ■■■■Xb)"^ and then carefully dividing numer- 
ical prefactors in appropriate proportions amongst rele- 
vant rate variables ij . Probability conservation demands 
this connection to the multinomial expansion, but it's the 
parceling out of terms that guarantees the complexity of 
the model. With this set of ODEs in hand, essentially 
providing a coarse-grained real-space renormalization- 
group [RSRG] prescription of the original Seceder Model, 
we follow the flow equations for the concentration vari- 
ables, characterizing all relevant fixed points. Within this 
broader mB space, the d = 1 Seceder Model exhibits its 
full richness. As an indication of the wealth of this geo- 
metric pattern formation, consider for the moment triplet 
selection dynamic, to = 3, where a self-similar hierarchy 
of multibranch fixed points emerges. We examine, to il- 
lustrate, the case B ^ 4, for which: 

Xl = xl + 3xi{x2 + -I- ^4) -I- 3x1X2X3 -\- 6X1X3X4 — Xl 
X2—X2+ 3x2{xl + ^3 -I- xl) + 3X2X32:4 — X2 

Because of branch symmetry about the central axis, the 
flow equations for the remaining variables are easily ob- 
tained via the interchange xi <-> x^ and X2 2:3 and, 
indeed, the globally stable fixed point (FP) must lie 
within this reduced subspace, with mirror variables iden- 
tified. For the case at hand, invoking the constraint 
X2 = 1/2 — 2:1 and demanding xi = leads to the cubic 
equation 7x5 — 6x4-1-5/4x1 = 0, yielding (xi, X2, X3, X4) = 
(r^, i, i, ^), as well as the less stable 2-branch solutions 
(^,0,0,^) and (0,^,^,0). Of course, if we simply nu- 
merically integrate the coupled ODEs and follow the tra- 
jectories from a randomly generated initial condition, we 
fiow with 100% probability to our unique superstable FP. 
The situation for i? = 5 is slightly different- insisting 
upon symmetry X5 — xi and X4 = X2, in addition to the 



normalization constraint = 1 — 2x2 — 2xi, we have the 
recurring 2-branch solutions with xi = and X2 = 1/2 
and vice versa, leaving us with two coupled bilinear equa- 
tions in the variables Xi and X2, represented graphically 
as a rotated, displaced ellipse and hyperbola within the 
unit square. There are two intersection points: one full- 
fledged superstable 5-branch solution, 
the other, an unstable lower-dimensional 4-branch so- 
lution (i, i, 0, 1;, i). Note that, in the latter instance, 

xi + X2 = 1/2 = Xi + X5, so that this unstable solution 
can, thanks to the gap X3 = 0, be understood, via coarse- 
graining, as literally self-similar to its two-branch cousin 
(^,0, ^). This sort of hierarchical connection manifests 
itself regularly whenever we uncover a lower-dimensional 
fixed point; i.e., vanishing Xi in the branch structure. 
The most compelling instance of this phenomenon ap- 
pears when we search for a superstable 8-branch FP. In 
fact, there is none. The ODE flows converge on a pe- 
cuhar 6-branch solution, (ii, i,0, j^, Jq,0, |, which 
is exactly self-similar to the strongly attractive 3-branch 
fixed point (|, i, |). Interestingly, this Cesser 8-branch so- 
lution is distinct from the straight out 6-branch, roughly 
(0.30,0.12,0.08,0.08,0.12,0.30), easily shown to be irra- 
tional, as is the 7-branch and all those beyond 8. The 
9, 10 and 12-branch FPs show no zeros, but such be- 
havior becomes increasingly rare. The 11-branch so- 
lution has two gaps, x^ — xs — 0, but with xi = 

xii « 0.271, X2 = xio ~ 0.089, X3 = xg « 0.044 and 
X5 = 2:7 « 0.047, ~ 0.100, can be coarse-grained to 
a broad 3-branch again, though in this case only ap- 
proximately, self-similar to our dominant FP 
likewise, the 15-branch, although 19 and 21-branch FPs 
show three gaps and a self-similarity to the 4-branch 
(A. L L I.) 
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An additional payoff of this RSRG treatment of the 
Seceder Model is an explanation of the relative stabil- 
ity of 2 and 3-branch solutions for triplet (m = 3) and 
higher multiplet (m > 4) selection groups. Recall Figure 
1, which made it clear that for large m, the dynamics of 
the d = 1 stochastic Seceder Model are controlled entirely 
by the strongly attractive fixed point of the 2-branch pat- 
tern. For triplet selection, however, a hierarchy of multi- 
branch solutions is manifest, which meant the frequent 
appearance of a 3-branches, and somewhat occasionally 
4 and 5-branches, but a complete absence of 2-branch 
behavior. The essential dichotomy can be understood 
graphically by following the flows for 3-branch dynamics 
in the deterministic Seceder Model, assuming triplet se- 
lection group TO = 3, illustrated in Figure 2a, where we 
show the [111] plane X1+X2 + X3 — 1. We find our super- 
stable fixed point, (|, i, |), within the unit triangle. We 
note, in particular, that the outmost 2-branch solution 
(^,0, ^) is unstable to small perturbations off the edge. 
In turn, the single branch FPs at the triangle vertices are 
entirely unstable. As m — > 4, however, this interior FP 
merges with that at the midpoint on the triangle's lower 
edge, reversing the flow and stabilizing the 2-branch dy- 
namic. From this vantage point, it is clear that quartet, 
rather than triplet, selection is the marginal fact 
quickly confirmed by a stability analysis of the (i, 0, i) 
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FIG. 2. a) RSRG flows in a:;ia;2i3— space for the d = 1 de- 
terministic Seceder Model. A superstable 3-branch FP exists 
within the equilateral triangle for m < 4 only, b) Middle 
branch growth rate, ±2, for different multiplet selection group 
sizes, m = 3 — 5, showing stability of the 2-branch solution 
(1/2,0,1/2) for quartet selection and greater. 



FP for arbitrary to. To linear order, we find X2 — 
[2(")(i)""^ - l]x2 + 0{xl), so we flow back to van- 
ishing X2 for 2m < 2'"~^; i.e., to > 4, since prefactor 
of the quadratic term is negative for the marginal value 
TOc — 4 : see Figure 2b, which shows the full behavior. 
In fact, for a continuously variable selection group of size 
TO — 4— £, the perturbatively stable 3-branch FP is located 
off the triangle's lower edge at X2 = r^(ln2 — j). 

As one considers the stability of higher {B > 3) multi- 
branch FPs, rric — 4 remains the upper critical size of 
the selection group at, and above which, mean-field deter- 
ministic pattern formation holds sway. For example, with 
quartet selection in the context of the 4-branch solution 
(i.e., TO = 4, B = 4), the superstable FP hes at the mid- 
point of the edge connecting the xi and Xi vertices of the 
unit tetrahedron- that is, {xi,X2,X3,X4) = (^,0,0,^). 
All initial starting points within the tetrahedron fiow 
outwards to the FP on this edge. While the other five 
edges of the tetrahedron are stable along their lengths, 
they're unstable in all other directions; the vertices, cor- 
responding to 1-branch solutions, are maximally unsta- 
ble. For triplet selection, the superstable 4-branch FP is, 
as mentioned earlier, (^, i, i, ^), which by contrast lies 
within the tetrahedron. There are also stable 3-branch 
FPs in this case located in the faces of the tetrahedron, 
where one of the Xi = 0, but these are unstable to pertur- 
bations toward interior of the tetrahedron. Vertex FPs 
are unstable to edge FPs, which in turn are unstable to 
face FPs, etc. Interestingly, these findings might suggest, 
at least initially, that the greatest stability is associated 
with the largest number of branches; i.e, {B — 1)— branch 
being unstable to _B— branch solutions, etc. However, 
the absence of superstable 8-branch FP, for starters, and 
the subsequent appearance of gaps for B = 11, 13, 17 in 
the spectrum indicate that any runaway tendency toward 
proliferation of branches from tip splitting will be cut off. 




FIG. 3. Despite greatly biased initial conditions, with two 
distinct localized groups, the d = 2 triplet Seceder model 
evolves the population toward an endpoint of three equidis- 
tant, separating clusters. Here, A''=2000, and the temporal 
snapshots, indicated by different colors, correspond to succes- 
sive times t — 2^, for i = 2 — 8; red— mauve, respectively. 




FIG. 4. Self-organized secession in d = 3, as well as higher 
dimensions, collapses to the plane defined by three divergent 
groups. Population size A'^=1000; generations t — 2^~^^ . 



Indeed, that is precisely the characteristic behavior of 
the triplet Seceder Model, where 3-branch dynamics are 
typically seen, with occasional 4, 5 or 6-branch runs. 

We should stress, in this regard, that the 3-branch so- 
lution is an extraordinarily robust feature of this model, 
becoming even more so in higher dimensions, where the 
genotype is specified by an d— component vector rather 
than a single real number, the case we've focussed on thus 
far. For example, in d = 2, an initially homogeneous, or 
highly polarized population for that matter- see Figure 3, 
will eventually segregate into three distinct groups head- 
ing off along the symmetry axes of a triangle, each cluster 
equidistant from the other two. In d = 3 dimensions, we 
might expect four groups, perhaps, localized at the cor- 
ners of an expanding tetrahedron, preserving the notion 
of equal distance. Interestingly, however, this does not 
happen at all. Again, we observe the formation of just 
three groups- note Figure 4; the effect is stronger still for 
d > 4. Apparently, asymptotic higher-dimensional seces- 
sion involves segregational collapse to a greatly reduced, 
two-dimensional, subspace- the hyperplane defined by 
three fuzzy points of an expanding equilateral triangle 
whose angular orientation may vary from one realization 
to the next, but whose essential geometry does not. In- 
terestingly, this dimensional reduction can be understood 
within the context of the deterministic model- one con- 
siders the stability of d -I- 1 equally separated groups in 
d dimensions to perturbations (ultimately, statistical in 
nature) that bring one group closer to the rest For 
d > 3, we find that the minority group is unstable and 
will go extinct, whereas for d < 3, the zero FP associ- 
ated with this vanishing group reverses stability, yielding 
three separatist clusters whose relative population sizes 
are set by the degree of symmetry breaking. 

In sum, we have revealed the mean-field limit of the 
multidimensional Seceder Model. For selection multiplet 
sizes m > = 4, the nonequilibrium dynamics produce 



a steady-state with two opposing groups, independent 
of d. In the extreme societal Seceder limit (m N), 
the noisy dynamics dies away, leaving two tightly knit 
groups. For m = 3, multiple groups are typical, with 
three the norm. Higher dimensional genotypes/strategies 
produce, surprisingly, no further fragmentation. Us- 
ing a coarse-grained RG prescription, which discretizes 
and renders deterministic the model, we analytically un- 
cover a self-similar hierarchy of multiple branch FPs in a 
gapped spectrum. Additional work, concerning the inter- 
mittent extinction dynamics of individual groups, early- 
time transient behaviors, as well as kinetic symmetry- 
breaking phenomena, will be reported elsewhere [ pO{ . 
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